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This  report,  the  fourth  in  a  series,  describes  the  techniques 
utilized  to  numerically  hlndcast  oceanic  wind  fields  and  presents  com¬ 
parisons  of  the  hindcast  wind  data  and  measured  wind  information.  The 
study  was  conducted  in  the  Hydraulics  Laboratory  under  the  direction  of 
Mr.  H.  B.  Simmons,  Chief  of  the  Hydraulics  Laboratory,  Dr.  R.  W.  Whalin, 
Chief  of  the  Wave  Dynamics  Division,  Mr.  C.  E.  Chatham,  Jr.,  Chief  of 
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Resio,  Dr.  C.  L.  Vincent,  and  Mr.  W.  D.  Corson.  Dr.  Resio  is  presently 
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OBJECTIVE  SPECIFICATION  OF  ATLANTIC  OCEAN 


WIND  FIELDS  FROM  HISTORICAL  DATA 


Introduction 


1.  The  single  most  critical  factor  in  calculating  an  oceanic  wave 
climate  by  means  of  numerical  models  is  the  accurate  estimation  of  wind 
fields  and  the  associated  momentum  exchange  from  the  atmosphere  into  the 
wave  field.  Report  1*  in  this  series  elaborated  the  procedure  for  ob¬ 
taining  representative  pressure  fields  from  historical  data.  This  re¬ 
port  discusses  the  methodology  employed  in  converting  these  pressure 
fields  into  wind  fields  at  the  surface  level.  It  is  important  to  point 
out  here  that  the  approach  used  attempts  to  simulate  the  physics  of 
synoptic-scale  atmospheric  motions  to  the  degree  permitted  by  our  under¬ 
standing  of  these  phenomena.  Consequently,  there  are  no  arbitrary  con¬ 
stants  that  can  be  adjusted  to  "fit"  the  data  in  individual  or  composite 
comparisons.  The  success  of  this  approach  as  indicated  in  the  many  com¬ 
parisons  to  be  presented  demonstrates  that  although  such  models  of  the 
physics  contain  various  assumptions  implicit  in  their  formulation,  the 
estimated  winds  should  be  representative  of  actual  wind  conditions  in 
any  oceanic  area.  The  possible  exception  to  this  would  come  in  tropical 
areas  where  subsynoptic-scale  phenomena  become  Important  and  the  geo- 
strophic  balance  assumption  becomes  a  poor  approximation. 

2.  The  type  of  wind  specification  used  here  plays  an  important 
role  in  the  reconstruction  of  past  individual  meteorological  and  ocean¬ 
ographical  events,  such  as  surface  waves  and  storm  surges  from  severe 
extratropical  cyclones.  Similarly,  in  a  probabilistic  treatment  of  ex¬ 
pected  patterns  of  surface  drift  and  associated  environmental  phenomena. 


*  W.  D.  Corson,  D.  T.  Reslo,  and  C.  L.  Vincent.  1980  (Jul).  "Wave 
Information  Study  for  U.  S.  Coastlines;  Surface  Pressure  Field  Recon¬ 
struction  for  Wave  Hindcasting  Purposes,"  Technical  Report  HL-80-11, 
Report  1,  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE, 
Vicksburg,  Miss. 
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the  dynamics  of  a  system  often  must  be  related  to  long-term  patterns  of 
winds,  waves,  currents,  and  tides.  It  is  frequently  not  sufficient  to 
substitute  a  mean  wind  field  into  such  a  climatological-scale  relation¬ 
ship.  Instead,  sequences  of  synoptic  events  with  their  wind  field 
structures  and  related  wind  speeds,  fetches,  and  durations  must  be  con¬ 
sidered.  In  viewing  this  problem  of  historical  prediction  for  a  wide 
range  of  potential  applications,  one  must  consider  not  only  contamina¬ 
tion  of  the  final  wind  field  due  to  modeling  errors  given  accurate  orig¬ 
inal  data  but  also  the  propagation  of  errors  in  the  original  data  through 
the  model  into  the  results.  Thus,  it  is  not  always  intuitively  clear 
whether  a  more  sophisticated  or  less  sophisticated  model  will  produce 
better  results  when  applied  to  actual  data,  containing  various  types 
and  magnitudes  of  errors. 

3.  Whereas  wind  fields  in  tropical  storms  have  been  represented 
in  many  studies  by  a  small  set  of  parameters,  the  complexity  and  diver¬ 
sity  of  shapes,  intensities  and  internal  structures  of  extratropical 
storms  have  precluded  the  use  of  parametric  wind  fields  for  these 
storms.  In  this  report,  tropical  storms  shall  be  excluded  and  attention 
restricted  to  the  representation  of  past  synoptic-scale  extratropical 
systems. 

4.  Four  primary  sources  of  information  were  included  in  the  as¬ 
similation  of  these  wind  fields: 

a.  Gridded  Northern  Hemisphere  pressure  fields  (Holl  and 
Mendenhall  1971)  (available  from  Navy  Fleet  Numerical 
Weather  Central,  Monterey,  Calif.). 

_b.  North  American  Historical  Weather  Map  Series  (available 
from  National  Climatic  Center,  Asheville,  N.  C.). 

£.  Surface  marine  observations  (TDF-11  available  from  Na¬ 
tional  Climatic  Center,  Asheville,  N.  C.). 

d_.  Airways  surface  observations  (TDF-14  available  from  Na¬ 
tional  Climatic  Center,  Asheville,  N.  C. ). 

Other  than  the  weather  maps,  all  of  this  information  is  available  on 
magnetic  tapes.  It  was  decided  to  restrict  attention  first  to  computer 
processible  forms  of  data  and  to  include  other  sources  of  data  if  it 
were  determined  that  these  could  lead  to  a  significant  improvement 


4 


in  the  accuracy  of  the  final  wind  field. 

5.  Basically,  two  independent  types  of  information  in  the  four 
data  sets  are  considered  here — pressures  and  winds.  The  pressure  data 
must  first  be  converted,  by  means  of  an  analysis  of  pressure  gradients, 
into  estimates  of  quasi-geostrophic  winds.  Then,  these  must  be  trans¬ 
formed  into  estimates  of  the  wind  vectors  at  some  reference  level  near 
the  surface.  Only  after  this  is  it  possible  to  combine  the  two  sets  of 
data,  using  some  blending  algorithm,  into  a  single  wind  field. 


Estimation  of  Quasi-Geos trophic  Winds 


6.  All  of  the  equations  used  in  this  report  for  the  estimation  of 
winds  at  the  geostrophic  level  can  be  found  in  any  basic  meteorology 
text  (for  example,  Hess  1959).  Consequently,  in  this  section,  the  math¬ 
ematics  will  be  omitted  with  the  understanding  that  the  equations  for 
simple  flows,  which  involve  a  balance  of  forces  (geostrophic  winds, 
gradient  winds) ,  baroclinic  effects  (thermal  winds) ,  and  response  to 
certain  pressure  field  changes  (isallobaric  winds),  can  be  found 
elsewhere. 

7.  Figure  1  shows,  on  a  polar  projection,  the  grid  system  on 
which  pressure  values  were  available.  A  preliminary  check  of  the  pres¬ 
sure  fields  in  this  data  set  and  those  in  the  weather  map  series  revealed 
some  large  discrepancies  between  the  two.  As  shown  in  Figure  2,  the 
gridded  data  appear  to  misrepresent  interiors  of  severe  extratropical 
lows.  A  reduction  of  intensity  of  these  systems  was  found  to  be  typical 
of  the  gridded  pressure  data  and  probably  represents  a  combination  of 
the  scale  of  the  grid  and  the  blending  procedure  used  to  smooth  the  pres¬ 
sure  data.  Actual  surface  observations  show  much  better  agreement  with 
the  weather  map  pressure  fields;  thus,  it  is  presumed  that  the  gridded 
data  are  in  error.  Given  that  this  is  the  case,  it  would  not  be  possi¬ 
ble  to  obtain  an  accurate  estimate  of  the  geostrophic-level  winds  from 
the  gridded  pressure  specification  without  some  modification. 

8.  In  the  western  Atlantic,  cyclogenesis,  rapid  storm  movement, 
and  complexity  of  storm  structure  make  it  difficult  to  find  a  simple 
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Figure  1.  Meteorological  International,  Incorporated  (Mil), 
63  x  63  grid  on  polar  projection  of  the  Northern  Hemisphere 
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method  for  modifying  the  gridded  data  into  pressure  fields  equivalent  to 
those  of  the  weather  maps.  For  the  area  shown  in  Figure  3,  weather  maps 


Figure  3.  Map  indicating  area  in  which  special 
attention  was  given  to  storms 

(available  at  6-hr  intervals)  were  reviewed  and  those  containing  at 
least  two  closed  contours  within  this  area  where  machine-digitized.  For 
storms  outside  of  this  area,  a  respecification  of  central  pressure  was 
used  to  amplify  the  pressure  gradients  in  the  storm  interior.  In  order 
to  maintain  resolution  of  steep  gradients,  these  respecified  areas  and 
digitized  areas  were  blended  into  an  interpolated  pressure  field  on  a 
grid  with  one-fourth  of  the  spacing  of  the  original  grid  (-30  NM) . 


Direct  Pressure  Field  Components 


9.  Calculations  of  the  surface  wind  fields  were  made  in  a  coordi¬ 
nate  system  that  consisted  of  great  circle  paths  for  quasi-east-west 
lines  and  orthogonals  to  the  great  circle  paths  for  quasi-north-south 
lines  (Figure  4).  The  reason  for  this  was  to  have  great  circle  paths 
for  wave  propagation  into  the  east  coast  of  the  United  States,  our  prin¬ 
cipal  area  of  concern  in  this  study.  The  basic  methods  could  be  applied 
to  other  coordinate  systems  after  appropriate  transformations. 

10.  The  first  approximation  to  the  surface  wind  field  in 


Figure  4.  Have  Information  Study's  (HIS)  spherical  orthogonal  grid 
on  a  Mercator  projection.  Heavy  line  is  numerical  boundary 
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non tropical  regions  is  the  geostrophic  balance  of  pressure  gradient  and 
Coriolis  force.  The  computer  technique  was  designed  to  search  for  the 
local  maximum  pressure  gradient  across  a  particular  grid  point  and  to 
base  the  geostrophic  wind  speed  and  direction  on  this  value.  There  are 
two  major  occasions  where  the  surface  wind  field  is  influenced  by  the 
pressure  field  by  components  of  a  nongeostrophic  nature.  In  tightly 
curved,  low-pressure  systems  the  effect  of  centrifugal  accelerations 
modifies  the  pressure  gradient-Coriolis  force  balance  to  produce  lower 
geostrophic  wind  speeds.  In  such  situations  a  gradient  wind  estimate 
was  substituted  for  the  geostrophic  wind  estimate.  Both  the  geostrophic 
and  gradient  wind  approximations  are  based  on  the  presumption  of 
balanced  forces.  Yet  there  are  many  times,  particularly  in  the  northern 
and  western  Atlantic,  when  the  time  rate  of  change  of  the  pressure  field 
is  substantial  and  negates  the  balanced  force  assumptions.  In  these 
situations  the  isallobaric  component  of  the  wind  field  is  estimated  and 
included.  A  final  consideration  in  the  estimation  of  winds  at  the  geo¬ 
strophic  level  is  the  effect  of  horizontal  temperature  gradients.  Such 
gradients  introduce  a  thermal  wind  component  into  geostrophic-level  winds. 
In  order  to  compensate  for  this  effect,  the  thermal  wind  equations  were 
integrated  up  to  a  nominal  level  of  1000  m  and  this  term  was  added  to 
the  resultant  of  the  previous  wind  components. 

Estimation  of  Winds  near  the  Surface 

11.  In  recent  years,  considerable  research  has  gone  into  comput¬ 
ing  interactions  between  the  free  atmosphere  and  the  planetary  boundary 
layer.  Since  the  scales  considered  in  numerical,  general  circulation 
models  cannot  resolve  fluxes  of  momentum,  heat,  and  moisture  from  the 
surface,  parametric  methods  are  used  in  these  models  to  couple  the  upper 
level  flow  to  the  earth's  surface  through  the  planetary  boundary  layer. 

The  most  successful  methods  of  parameterizing  fluxes  in  the  boundary 
layer  have  evolved  by  means  of  matching  solutions  from  similarity 
theories  for  an  outer  layer  (Ekman  region)  and  inner  layer  ("constant 
stress"  region) . 


12.  In  a  diabatic,  baroclinic  atmosphere  Hess  (1973)  has  shown 
that  the  flow  in  the  boundary  layer  follows  a  nondimens ional  geostrophic 
departure  function  of  the  form 


(1) 
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where  k  is  the  von  Karman  constant,  a  solution  for  the  winds  in  the 
boundary  layer  can  be  obtained  in  the  form  (Clarke  and  Hess  1974) 


la  Ro  ■  A  (a,  Sx>  Sy)  - 


In  — 


k2|  V 


-|l/2 


-  B 


(“•  V  Sy) 


(2) 
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sin  0 


k  V 


where 


surface  Rossby  number 


(3) 


zq  =  surface  roughness  length 

■+ 

V  =  surface  geostrophic  wind 
8o 

0  *  angle  between  the  surface  geostrophic  wind 
and  the  surface  stress 

A,B  =  two  nondimens ional  functions 


13.  Recent  evidence  (Kitaigorodskii  and  Volkov  1965,  Smith  and 
Banke  1975,  Garratt  1977)  suggests  that  zq  is  not  an  external  param¬ 
eter  for  a  water  surface,  but  rather  is  an  implicit  function  of  momentum 
fluxes  near  the  surface.  Cardone*  has  revised  his  functional  form  of 
zq  (Cardone  1969)  to  be  in  agreement  with  the  findings  of  Garratt  (1977) 

z  *  0.1525  uT1  +  0.0144  u 2  -  0.00371  (4) 

O  "  ” 


(in  cgs  units).  Equations  2,  3,  and  4  along  with  an  equation  for  the 
mean  velocity  profile  in  the  "constant  stress"  region, 

_ -Ui  -  [*-  (r)  -  ♦(!)] _ 

*  Personal  communication,  1977. 
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where  <j>  is  a  universal  similarity  function,  provide  a  system  of  equa¬ 
tions  that  can  be  solved  iteratively  for  u^  ,  z0  »  an<*  ®  •  The 
forms  of  \p  used  in  Equation  5  are  the  KEYPS  formula  (Lumley  and 
Panofsky  1964)  for  unstable  stratification  and  a  log-linear  function 
for  stable  stratification.  Details  of  how  typically  recorded  air  tem¬ 
perature  and  sea  temperature  can  be  incorporated  into  these  equations 
can  be  found  in  Cardone  (1969). 

14.  The  functions  A  and  B  in  Equations  2  and  3  contain  a  de¬ 
pendence  on  both  stability  and  geostrophic  shear.  Arya  and  Wyngaard 
(1975)  demonstrated  that  these  two  factors  can  be  considered  indepen¬ 
dently  and  hence  A  and  B  can  be  written  as 


A  "  As  +  S 
B  =  Bs  +  Bb 


(6) 


where  the  subscripts  s  and  b  refer  to  effects  of  stability  and  baro- 
clinicity,  respectively.  For  a  linearly  decreasing  shear  to  the  top  of 
a  boundary  layer  scaled  by  u^/f  ,  the  Afa  and  Bb  terms  integrated 
over  the  boundary  layer  thickness  become 

A^  =  k  cos  0 

isi  <7> 

Bb  *  k  -Jy1-  sin  0 

where 

3  =  dimensionless  thermal  vector 

6  ■  angle  between  the  geostrophic  wind  and  this  vector 
The  operator  3  is  used  here  to  define  the  horizontal  temperature 
gradient  at  each  point.  The  forms  of  the  Ag  and  Bg  terms  are  param¬ 
eterized  in  terms  of  the  dimensionless  stability  parameter  p  , 
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(8b) 


Figure  5  shows  a  comparison  of  these  parameterlzations  to  field  data 
from  Clarke  and  Hess  (1974) .  As  seen  here  for  p  >  0  ,  the  A  and  B 
functions  appear  to  approach  asymptotic  limits  of  about  +15  in  the 
field  data;  consequently,  these  limits  were  used  in  the  PBL  computer 
subroutine  to  replace  calculated  values  greater  than  +15.  In  this 
range  the  solution  for  the  surface  stress  became  insensitive  to  small 


Figure  5.  Comparison  of  parameterlzations  to  field 
data  (from  Clarke  and  Hess  1974) 


variations  in  A  and  B  so  the  exact  values  become  somewhat  unimpor- 
s  s 

tant.  This  assumption  ensures  that  no  unreasonable  magnitudes  of  p  come 
into  the  PBL  routine  and  thus  avoids  potential  problems  with  the  numeri¬ 
cal  stability  of  the  computations. 

15.  It  should  be  noted  that  in  this  report  the  scale  height  used 
in  the  parameterization  of  the  boundary  layer  is  u^/f  and  not  the 
height  of  the  inversion  layer.  Although  this  might  lead  to  serious  mis¬ 
representations  in  boundary  layers  over  land,  over  an  open-ocean  area  it 
appears  justifiable  to  retain  a  u^/f  scaling,  due  to  much  lower  rates 
of  surface  heating  and  cooling.  As  with  most  assumptions  of  this  nature, 
only  results  of  extensive  comparisons  will  demonstrate  whether  or  not  it 
is  reasonable. 

16.  Given  a  specification  of  stability  in  terms  of  a  temperature 
deviation  between  the  water  surface  and  the  air  at  some  reference  level, 
the  stability  length  can  be  estimated;  and  given  a  value  for  the  horizon¬ 
tal  temperature  gradient  in  the  surface  layer,  the  thermal  wind  effects 
can  be  calculated.  It  should  be  immediately  obvious  that  the  quality 

of  such  estimates  can  be  criticized.  In  fact,  a  primary  concern  in  in¬ 
troducing  a  dependence  on  stability  and  baroclinicity  into  the  wind  es¬ 
timation  methodology  is  the  possibility  of  large  errors  in  the  input 
thermal  fields.  However,  rather  than  assume  a  priori  that  the  informa¬ 
tion  content  of  these  estimates  will  negatively  affect  the  accuracy  of 
the  wind  calculations,  it  was  decided  to  proceed  on  the  hypothesis  that 
the  overall  error  in  wind  estimates  would  be  reduced  by  inclusion  of 
some  measure  of  the  horizontal  and  vertical  thermal  gradients  in  the 
planetary  boundary  layer. 

17.  In  Appendix  A  a  derivation  of  a  comparison  of  inherent  errors 
in  wind  fields,  with  and  without  consideration  of  thermal  factors,  is 
given.  As  indicated  there,  the  less  sophisticated  approach  does  not 
guarantee  better  results  on  the  basis  of  a  lack  of  sensitivity  to  uncer¬ 
tain  input  parameters.  Instead,  the  lack  of  consideration  of  persistent 
factors,  such  as  the  strong  temperature  gradients  along  the  U.  S.  Atlan¬ 
tic  coast,  can  introduce  large  regional  biases  in  wind  estimates.  On 
the  other  hand,  it  is  also  important  to  avoid  the  introduction  of  large 


15 


local  errors  Into  objective  schemes  through  the  use  of  overly  sensitive 
measures  of  thermal  gradients.  For  example,  if  there  is  a  sparse  network 
of  ships  reports  of  surface  temperatures  and  air-sea  temperature  differ¬ 
ences,  a  single  erroneous  ship  report  can  introduce  significant  errors 
into  an  estimate  of  thermal  gradients  over  a  relatively  large  area. 

Also,  there  will  be  times  in  some  regions  where  no  reliable  information 
exists  for  a  prolonged  period  of  time.  Since  the  hindcast  effort  for 
which  these  winds  were  needed  extends  over  a  20-year  interval  (1956 
through  1975) ,  the  effort  required  to  remedy  such  situations  would  have 
been  almost  insurmountable.  Thus,  a  different  approach  was  needed. 

18.  After  some  analyses,  it  was  determined  that  stable  estimates 
of  local  surface  temperature  and  air-sea  temperature  difference  could 

be  obtained  by  stratifying  their  variation  with  respect  to  time  into  two 
scales,  seasonal  and  synoptic.  Time  of  year,  t*  ,  sufficed  to  give  a 
good  measure  of  the  former  scale  of  variation;  and  the  direction  of  the 
geostrophic  wind,  a  ,  provided  a  reasonable  measure  of  the  flux  of  air 
from  different  source  regions  on  a  synoptic-scale.  Hence,  the  air  tem¬ 
perature,  T  ,  and  air-sea  temperature  difference,  T  ,  could  be 
characterized  by  two  parameters.  To  obtain  a  reliable  estimate  of  T^ 
and  T  _  ,  the  total  sample  of  ships  observations  for  four  5-year 

fl“8 

periods  was  used  to  define  these  variables  for  each  month  of  the  year  for 
each  10-deg  wind  direction  category.  For  each  oceanic  grid  point  shown 
in  Figure  4,  these  two-dimensional  matrices  were  calculated.  Each  ther¬ 
mal  field  constructed  in  this  manner  (as  a  function  of  t*  and  a)  was 
smoothed  after  the  initial  estimate  to  further  reduce  the  possibility  of 
extreme  gradients  over  one  or  two  grid  points. 

19.  After  the  geostrophic-level  winds  were  reduced  to  near-surface 
winds.  Independent  observations  from  ships  were  blended  into  the  wind 
fields.  The  procedure  used  was  simply  to  average  the  observed  wind 
speed  with  the  nearest  grid  point  wind  speed.  Although  more  sophisti¬ 
cated  methods  were  considered,  the  effects  due  to  differences  in  storm 
position  seemed  to  offset  any  potential  gain  inherent  in  a  more  elabo¬ 
rate  methodology.  Wind  directions  from  ships  were  not  blended  into  the 
wind  fields  since  a  180-deg  direction  error  can  be  quite  common  due  to 
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the  manner  in  which  the  values  are  coded.  Over  the  20-year  period  for 
which  calculations  were  made,  approximately  7.4  million  ships  observa¬ 
tions  were  averaged  into  the  final  wind  fields. 

Comparison  of  Results 

20.  In  order  to  investigate  whether  or  not  the  wind  fields  being 
calculated  provide  a  reasonable  representation  of  actual  conditions,  ex¬ 
tensive  comparisons  between  predicted  winds  and  independent  observations 
have  been  made.  Three  primary  sources  of  data  were  involved  in  these 
comparisons:  recorded  data  from  NDBO  (NOAA  (National  Oceanic  and  Atmos¬ 
pheric  Administration)  Data  Buoy  Office)  buoys.  Sable  Island  anemometer, 
and  ships  observations  (before  th**  information  from  these  observations 
was  blended  into  the  wind  fields).  Along  the  Atlantic  coast,  hourly 
measurements  of  wind  speed  and  direction  at  the  10-m  elevation  from  the 
buoys  afford  an  excellent  opportunity  to  check  the  details  of  winds 
through  time  at  a  specific  point.  To  establish  the  validity  of  the  wind 
estimation  procedure  over  the  entire  Atlantic  Ocean,  analyses  must  be 
made  on  the  basis  of  comparisons  within  different  portions  of  the  Atlan¬ 
tic  Ocean.  The  only  data  presently  available  that  can  provide  good 
coverage  over  the  entire  Atlantic  Ocean  are  ships  observations.  Unfor¬ 
tunately,  most  of  the  observations  made  in  the  Atlantic  Ocean  are  visual 
estimates  of  wind  speeds  and  directions  rather  than  instrument  measure¬ 
ments;  therefore,  any  comparisons  made  to  them  must  consider  the  effects 
of  biases  and  random  errors  in  the  visual  observations. 

21.  As  previously  stated,  the  available  observed  data  sets  are 
different  in  nature,  i.e.,  the  ships  data  are  for  relatively  long  periods 
of  time  and  cover  the  Atlantic  Ocean  (nonspecific),  while  buoy  or  plat¬ 
form  data  are  more  specific  in  time  and  space  (specific).  The  nonspeci¬ 
fic  ships  data  will  be  compared  in  this  section  and  comparisons  to 
specific  information  will  be  presented  in  Appendix  B.  Also,  the  compari¬ 
sons  to  ships  wind  speed  will  be  discussed  first,  and  then  the  compari¬ 
sons  of  wind  direction  will  be  presented. 

22.  Before  an  estimate  of  the  random  error  can  be  obtained,  bias 
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present  in  the  ships  observations  must  be  removed  to  whatever  extent 
possible.  Cardone  (1969)  summarizes  previous  investigations  of  bias  in 
winds  reported  by  ships  by  the  following  formula  (in  knots) : 

w  =  2.16  w7/9  (9) 

r 


where 

w  =  corrected  wind  speed 
w^  =  reported  wind  speed 

Unless  otherwise  stated,  this  relationship  was  used  to  transform  ships 
wind  speeds  prior  to  comparison  to  predicted  wind  speeds. 

23.  The  magnitude  of  random  errors  in  ships  observations  of  wind 
speed  can  be  estimated  from  comparisons  to  some  known  standard  such  as 
an  NDBO  buoy.  Figure  6  shows  the  results  of  such  a  comparison.  As  is 
typical  of  this  type  of  environmental  data,  the  random  error  is  not 
homoscedastic,  but  rather  has  a  dependence  on  the  magnitude  of  the 
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variate.  Consequently,  if  a  comparison  between  our  predicted  wind  speeds 

2 

and  ships  observations  is  made,  the  apparent  error  variance,  a  must 

be  stratified  by  velocity  class  in  order  to  be  able  to  estimate  its  two 

2 

components,  error  due  to  incorrect  predictions,  o  ,  and  error  due  to 

2  p 

inaccurate  observations,  oq  .  The  estimate  of  error  in  predicted 
winds  obtained  in  this  fashion  can  be  written  as 


f  2  2  ~\ 

ap(u)  =  Ua(u)  -  oq(u)  (10) 

2 

where  o  ,  determined  independently  from  comparisons  of  ships  winds  to 

°  2  —  — 
buoys,  is  approximated  by  o  *  1.5  +  0.18  u  ,  where  u  is  the  sample 

mean  of  the  wind  speeds. 

24.  From  comparisons  to  buoys  and  using  Equation  10  for  compari¬ 
sons  to  ships  observations  it  is  seen  that  there  is  an  apparent  depen¬ 
dence  on  wind  speed  of  the  standard  deviation  of  the  error  distribution 
which  is  approximately  by  the  linear  relation 


op  =  1.7  +  0.11  u  (11) 

where 

u  *  average  recorded  wind  speed  in  the  comparisons 
o  =  standard  deviation  of  the  errors 
This  error  compares  well  with  reported  error  values  in  Overland  and 
Gemmill  (1977).  In  that  study  a  limited  area,  fine-mesh  model  was  used 
to  obtain  pressures  and  wind  speeds  at  the  geostrophic  level  over  the 
New  York  Bight.  The  best  results  for  surface  winds  were  obtained  with 
a  PBL  model  (Cardone  1969)  approach,  producing  o  values  of  2.92  m/sec, 
3.49  m/sec,  and  2.93  m/sec  for  u  <  5,  5  <  u  <  10,  and  u  >  12.5  m/sec 
wind  speed  classes,  respectively. 

25.  Figures  7,  8,  and  9  present  three  separate  comparisons  of 
ships  and  predicted  wind  speeds.  The  distribution  of  deviations  is  dis¬ 
played  in  Figure  7.  Figure  8  shows  contoured  densities  for  a  cross  plot 
of  predicted  wind  speeds  and  observed  ships  winds  before  correction 
through  Equation  9  for  all  cases  in  1974.  Comparison  parameters  from 
the  January,  April,  July,  and  October  months  were  calculated  from  the 
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7.  Deviation  of  computed  wind  speeds  (grid)  from  measured 
wind  speeds  (ships)  for  one  year 


Figure  8.  Joint  distribution  of  paired  observations  of  observed 
(ships)  wind  speeds  and  predicted  (WIS)  wind  speeds 
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9.  Bias  and  root-mean-square  error  (RMSE)  of  wind  speed  for  the 
indicated.  Information  was  computed  for  10  years  of  each  month 
(some  dots  overlay  each  other) 


20  years  of  information  and  are  presented  in  Figure  9.* 

26.  Since  the  results  show  the  same  type  of  observer  bias  as  given 
independently  by  Equation  9,  they  indicate  that  there  is  little  apparent 
bias  in  the  wind  speed  estimates  across  the  Atlantic  Ocean. 

27.  As  shown  in  Figure  10,  there  does  not  appear  to  be  signifi¬ 
cant  bias  in  predicted  wind  directions  when  compared  with  ships  observa¬ 
tions.  Thus,  the  wind  field  estimation  methodology  outlined  in  this 
paper  seems  to  provide  a  good  synoptic-scale  description  of  actual  wind 
fields  over  the  entire  Atlantic  Ocean. 

28.  It  should  be  noted  that  in  our  comparisons  the  random  error 
in  wind  directions  decreases  markedly  with  increasing  wind  velocity. 

Upon  inspection,  it  appears  that  much  of  the  deviation  between  observed 
and  predicted  wind  speeds  and  directions  is  due  to  differences  in  storm 
locations.  Time-series  plots  of  simulated  winds  often  exhibited  lags 
or  leads  of  6  hours  or  so  compared  with  nearby  measurements.  Since 
winds  tend  to  be  less  well  organized  under  low  wind  speed  conditions, 
and  hence  more  highly  variable,  it  is  expected  that  such  conditions 
will  produce  larger  errors  in  wind  direction.  An  additional  factor  to 
consider  in  these  comparisons  is  that  nonins trumented  estimates  of 
wind  direction  appeared  to  be  highly  dependent  of  wind  speed. 

29.  It  was  mentioned  near  the  beginning  of  this  report  that  we 
would  reconsider  our  decision  to  ini  lude  only  large  data  sets  in 
machine-processible  form,  based  on  our  final  result.  Our  primary  in¬ 
terest  in  this  study  is  a  region  along  the  east  coast  of  the  United 
States,  where  a  high  density  of  ships  observations  prevailed  during  the 
20-year  interval  of  interest.  After  reviewing  the  results  obtained  in 
this  methodology  and  considering  the  large  number  of  data  already  in¬ 
cluded,  it  is  doubtful  that  any  small  set  of  additional  data  will 
strongly  influence  the  estimated  wind  fields.  Consequently,  the  deci¬ 
sion  was  made  not  to  process  additional  data  for  inclusion  into  the  wind 
fields  as  a  part  of  this  study.  This  is  not  intended  to  imply  that  some 


*  The  years,  from  which  the  months  were  chosen,  were  randomly  selected. 
The  parameters  of  each  averaged  month  are  intended  to  approximate  sea¬ 
sonal  information. 
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Figure  10.  Bias  and  RMSE  of  predicted  wind  direction.  Information  was 

selected  as  in  Figure  9 


local  Improvement  In  these  wind  field  estimates  may  not  be  obtainable 
by  blending  in  some  specialized  data  source  not  included  here. 


Summary  and  Conclusions 

30.  An  objective  method  for  obtaining  wind  fields  has  been  devel¬ 
oped  and  20  years  of  wind  fields  covering  most  of  the  North  Atlantic 
Ocean  have  been  calculated.  This  method  uses  historical  information 

on  pressure  fields  and  horizontal  thermal  gradients  to  obtain  an  esti¬ 
mate  of  the  wind  vector  at  the  geostrophic  level.  A  parametric  repre¬ 
sentation  of  the  planetary  boundary  layer,  based  on  similarity  theory, 
is  used  to  reduce  this  wind  to  a  surface  wind  vector.  Input  parameters 
in  terms  of  atmospheric  stability  and  baroclinicity  are  developed  from 
a  statistical  analysis  of  their  dependence  on  the  wind  at  the  geostrophic 
level.  Comparison  of  predicted  winds  to  observed  winds  indicates  that 
the  methodology  used  here  produces  relatively  unbiased  winds  with  a 
scalar  RMSE  slightly  dependent  on  wind  speed  (o^  »  1.7  +  0.11  u  m/sec 
where  u  is  the  mean  wind  speed  of  the  sample) . 

31.  The  hindcast  wind  fields  have  been  used  in  a  wave  hindcast 
program  to  generate  wave  data  for  the  North  Atlantic  Ocean  (Corson  et  al. 

1980)  .  Comparisons  of  the  hindcast  wave  data  and  measured  wave  data 
indicate  good  agreement  between  the  two  data  sets  (Corson  and  Resio 

1981) .  Thus,  it  is  apparent  that  the  methods  employed  to  generate  the 
20  years  of  wind  fields  were  appropriate. 
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Appendix  A:  Comparison  of  Expected  Errors  in  Wind  Estimation 
With  and  Without  Consideration  of  Stability  and 
Baroclinicity  in  the  Boundary  Layer 

1.  Neglecting  subsynoptic-scale  perturbations  in  the  wind  field, 
recent  experimental  evidence  and  theoretical  treatments  indicate  that 
the  wind  near  the  surface  should  be  primarily  a  function  of  five  factors, 
the  geostrophic  wind  vector  (£),  the  thermal  wind  vector  (V^T) ,  vertical 
temperature  gradients  near  the  surface  (V^T) ,  the  Coriolis  acceleration 
(f ) ,  and  the  roughness  of  the  surface  (zq)  (Cardone  1969*) .  A  "true" 
functional  relation  for  the  surface  wind  vector  thus  could  be  put  in  the 
form 

w  -  w  (6,  V^t,  Vl,  f,  zQ)  (Al) 

In  looking  for  an  optimal  wind  estimator,  we  are  generally  looking  to 
minimize  two  error  terms,  bias  and  random  error, 

b  -  E  [:  -  c*2> 

a2  *  E  [(w  -  t/]  (A3) 

where 

2 

b,a  »  expected  bias  and  error  variance,  respectively 

E  [  1  *  expectation  operator 

w  *  estimated  wind  vector  near  the  surface 
e 

Given  two  estimators,  one  which  considers  only  a  dependence  on  Z  and 
f  and  a  second  which  considers  all  five  primary  factors,  i.e. 

»i  -  wx  (<$,  f)  (A4) 

w2  -  w2  (Z,  rf,  f,  zo)  (A5) 

2 

then  b  and  a  can  be  defined  in  terms  of  integrals  over  a 
*  See  References  at  end  of  main  text. 


Al 


multidimensional  probability  space.  If  the  relation  for  w^  is  con¬ 
sidered  to  be  a  "true"  relation  for  neutral  conditions,  we  have  for  the 
case  ^G,  :G^=  0 

bl  =  f  f  <=  -  P  (tG.  V?.  *„>  <«£f)  d  <V>  d*0  <*« 

°i  »  (”  -  “!>2  p  <v  v^.  J0>  d‘cd  d  <V>  dzo  (A7> 

where  is  a  random  vector  describing  the  error  between  a  specified 

value  of  <5  and  the  true  value  of  G  .  If  we  allow  a  similar  set  of  ran¬ 
dom  error  vectors  to  characterize  the  errors  of  four  of  the  input  param¬ 


eters,  with  cv  ^o,  ez  0  , 

and  assume  that  the  underlying  relation  for  is  approximately  correct, 

then  we  obtain  for  b2 

b2  m  J  J  ‘  V  p  ( V  \r  \r  \)  dtGd\TdtvzTd*zo  (A8) 


and  for  a„ 


J...J  (w  -  w2)2  p  ^eG,  e 


V.T’  GV  T’  Ez  )  deGdEV  TdeV  Tdez 
h  z  o/  h  z  o 


The  important  distinction  between  the  integrals  in  A6  and  A7  and  A8 

and  A9  is  that  the  former  are  defined  over  ranges  in  input  variables 

(dV.  f ,  dV  i,  dz  )  as  well  as  over  a  range  in  the  random  vector  (de_) , 
n  z  o  g 

whereas  the  latter  are  defined  over  the  ranges  only  of  the  random  vectors 

(deG,  de^  T,  de7  T,  dez  J  .  Thus,  the  estimates,  which  are  typical 

of  winds  predicted  by  regression  method  (i.e.  |w|  ■  A | G | ) ,  where  A  is  a 
constant  for  a  specific  latitude,  can  contain  large  biases  when  applied 


A2 


in  regions  where  the  means  V,  T  and  V  T  are  not  equal  to  zero  or 

il  z 

where  z  is  different  from  the  z  for  which  the  statistical  method 
o  o 

was  calibrated.  For  example,  over  the  entire  Atlantic  Ocean,  the  average 
stability  is  near  neutral  and  the  average  baroclinicity  is  small;  how¬ 
ever,  near  the  east  coast  of  North  America  strongly  baroclinic,  unstable 

conditions  persist  during  the  winter. 

2 

2.  Similarly,  the  term  contains  contributions  from  the  ne¬ 
glect  of  the  ,  V  ^  ,  zq  terms.  If  the  dispersions  of  the  distri¬ 
butions  of  these  variables  ( E  ~  ^h^^l  *  etc-^  are  larger  than 


random  errors  in  the  specification  of 


w„ 


(E[vl  • etc) 


then  we 


should  have  o?  >  .  Since  the  variations  in  V,  "f  and  V  7  include 

12  h  z 

contributions  from  synoptic  scales  up  through  seasonal  scales  of  atmo¬ 
spheric  circulation,  they  can  indeed  be  quite  large  at  a  given  site. 

3.  On  the  other  hand,  if  the  dependence  of  on  the  input  vari¬ 

ables  is  almost  linear  or  if  the  dependence  is  symmetric  with  respect  to 
a  perturbation  around  a  "true"  value,  i.e. 


(-D 


v +  (-i)k%t-  + 

n 


(-DkV?’ f .  zo +  <-1)k%  1 

z  oj 

b 2  should  be  quite  small.  This  condition  is  very  nearly  met  for  the 
boundary  layer  operators  for  typical  e-values  .  Since  in  the  boundary 
layer  parameterization  here,  zq  is  an  implicit  function  of  the  surface 
stress  its  error  properties  are  included  in  the  effects  of  errors 

in  specifications  of  £  .  This  brings  out  an  interesting  point  since  a 
jwj  -  X|S|  relation  does  not  contain  this  implicit  dependence  on  zq  ; 
hence,  the  assumption  made  earlier  that  such  a  form  could  represent  a 
"true"  relation  for  neutral  conditions  is  not  strictly  valid. 

4.  Although  the  method  used  here  to  compare  expected  bias  and 


random  errors  seems  laborious,  it  does  provide  a  means  of  addressing  a 
somewhat  complex  question.  The  use  of  overly  simple  methodologies  is 
often  appealing  in  terms  of  convenience;  however  as  indicated  here, 
methods  which  better  represent  boundary  layer  effects  should  consistently 
provide  a  less-biased  estimate  of  the  surface  wind.  Also,  under  most 
circumstances,  the  PBL  approach,  including  thermal  effects  and  an  im¬ 
plicit  dependence  of  zq  on  the  surface  stress,  should  produce  less 
random  error  than  the  neglect  of  these  factors. 


Appendix  B:  Site-  and  Time-Specific  Comparisons 


1.  In  a  previous  section  of  this  report  extensive  comparisons  to 
ships  observations  were  presented.  These  comparisons  showed  the  pre¬ 
dicted  winds  to  be  valid  over  the  entire  Atlantic  Ocean.  In  this  sec¬ 
tion,  site-  and  time-specific  comparisons  are  employed  to  further  inves¬ 
tigate  characteristics  of  the  wind  hindcast  procedure. 

2.  Reliable,  recorded  wind  data  were  available  from  three  sites 
when  these  comparisons  were  performed*:  Sable  Island,  NDBO  buoy  EB01, 
and  NDBO  buoy  EB13  (Figure  Bl).  Data  from  each  site  were  compared  with 


Figure  Bl.  Location  map  for  available  recorded  data  site 


*  All  comparisons  have  been  adjusted  to  comparable  anemometer  heights. 


Bl 


data  from  the  nearest  Wave  Information  Study  (WIS)  grid  point. 

3.  The  wind  comparisons  at  Sable  Island  were  performed  after 
first  transforming  the  Sable  Island  winds  to  the  hypothetical  19.5-m 
level  of  the  WIS  winds.  Figures  B2-B6  indicate  that  the  wind  hindc-i^L 
data  are  in  phase  with  the  measured  data  at  Sable  Island  and  the  hind- 
cast  procedure  produces  winds  of  similar  magnitudes  as  the  measured  data. 
Additionally,  Figures  B7  and  B8  show  no  apparent  bias  in  the  predicted 
winds  at  Sable  Island. 

4.  To  better  describe  how  well  the  hindcast  procedure  performed 
for  high  wind  conditions,  computed  and  measured  cases  were  compared  for 
each  event  in  which  wind  speeds  over  30  knots  were  recorded  at  Sable 
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DAYS  IN  JANUARY,  t«72 


Figure  B2.  Comparison  of  hindcast  (WIS)  and  measured  (Sable  Island) 
wind  speed  at  Sable  Island  for  January  1972 
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Figure  B3.  Comparison  of  hindcast  (WIS)  and  mea¬ 
sured  (Sable  Island)  wind  speed  at  Sable  Island 
for  February  1972 
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DAYS  IN  MARCH,  1072 


Figure  B4.  Comparison  of  hindcast  (WIS)  and  measured  (Sable  Island) 
wind  speed  at  Sable  Island  for  March  1972 


B3 


DAYS  IN  NOVEMBER.  1972 


Figure  B5.  Comparison  of  hindcast  (WIS)  and 
measured  (Sable  Island)  wind  speed  at  Sable 
Island  for  November  1972 


Figure  B6.  Comparison  of  hindcast  (WIS)  and 
measured  (Sable  Island)  wind  speed  at  Sable 
Island  for  December  1972 
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Figure  B7.  Cross  plot  of  hindcast  (WIS)  and  measured  (Sable 
Island)  wind  speeds  at  Sable  Island  for  January  1972 
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WIND  SPEED,  KNOTS  (WIS) 


Figure  B8.  Cross  plot  of  hindcast  (WIS)  and  measured  (Sable 
Island)  at  Sable  Island  for  March  1972 


Island  (Figure  B9).  As  shown  in  Figure  B9,  the  computed  high  winds  are 
very  similar  to  the  measured  wind  speeds,  and  the  means  of  both  data  sets 
are  nearly  identical. 

5.  In  Appendix  A,  a  somewhat  heuristic  argument  was  presented 

in  favor  of  using  a  complete  simulation  of  the  planetary  boundary  layer 
in  the  estimation  of  surface  winds  from  geostrophic-level  winds.  A  past 
alternative  to  this  has  been  to  use  simple  algorithms  to  represent  such 
wind  relations.  Typically,  this  approach  is  justified  on  the  grounds 
that  a  more  thorough  treatment  is  not  merited  due  to  inherent  errors  in 
the  input  specifications.  As  shown  in  Appendix  A,  such  an  argument  is 
not  very  persuasive,  since  any  additional  bias  introduced  by  an  incor¬ 
rect  wind  model  will  add  to  any  existing  random  errors  and  biases  in  the 
wind  parameters  specified  at  the  top  of  the  boundary  layer.  The  ultimate 
goal  must  be  to  achieve  maximum  accuracy  in  the  specification  of  the 
geostrophic-level  winds  and  the  derived  surface-level  winds.  Figure  BIO 
presents  a  comparison  of  the  characteristic  behavior  of  a  similarly 
based  planetary  boundary  layer  model  to  a  simple  linear  algorithm  that 
is  still  used  in  different  engineering  and  scientific  applications 
today.  The  factor  "R"  referenced  in  this  figure  is  somewhat  a  "tuning 
coefficient"  and  usually  the  value  selected  is  referenced  as  being 
slightly  greater  than  1.  As  shown  in  Figure  BIO,  one  would  expect  that 
a  linear  relation  will  significantly  overpredict  extreme  wind  conditions. 
To  test  this  hypothesis,  winds  off  magnetic  tapes  from  Fleet  Numerical 
Weather  Central  (FNWC)*  were  compared  with  the  same  measured  winds  at 
Sable  Island  as  were  computed  winds  from  the  present  study.  The  hind- 
cast  winds  in  this  comparison  were  estimated  using  a  simple  algorithm 
rather  than  a  complete  boundary  layer  model.  As  shown  in  Figure  Bll, 
there  is  a  significant,  consistent  overprediction  at  high  wind  speeds. 

6.  Comparisons  to  the  NDBO  buoys  (EB01,  EB13)  contain  wind  speed 
and  direction  parameters.  Some  of  the  buoy  data,  although  plotted  as 
valid  data,  probably  represent  equipment  problems,  especially  some  of 

*  The  methodology  employed  to  compute  the  winds  for  FNWC  was  developed 
by  Meteorological  International,  Incorporated  (Mil)  and  is  discussed 
in  Holl  and  Mendenhall  1971. 
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SABLE  ISLAND  WIND  SPEED.  KNOTS 


Figure  B9.  Cross  plot  of  time-paired  hindcast  (WIS)  and  measured 
(Sable  Island)  wind  speeds  from  January  and  March  1972  for  obser¬ 
vations  at  Sable  Island  exceeding  30  knots 
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Figure  BIO.  Surface  wind,  Ug  ,  as  a  function  of  geostrophic  wind, 

Ug  ,  PBL  theory  and  from  a  simple  linear  algorithm  with  various 

tuning  coefficients  (R) 


SABLE  ISLAND  WIND  SPEED.  KNOTS 


Figure  Bll.  Cross  plot  of  hindcast  (Mil)  and  measured  (Sable  Island) 
wind  speeds  from  January  and  March  of  1972  for  observations  at  Sable 

Island  exceeding  30  knots 


BIO 


the  zero-wind  conditions.  In  spite  of  this,  the  time  series  plots  of 
predicted*  (WIS)  and  measured  (buoy  EB01  or  EB13)  wind  speeds  show  the 
predicted  winds  to  be  very  similar  to  the  measured  winds.  Likewise,  the 
time  series  plots  of  direction  differences  indicated  expected  variance 
with  most  of  the  larger  differences  occurring  at  low  wind  speeds 
(Figures  B12-B20) . 

7.  Figure  B21  presents  comparison  statistics  of  measured  (buoy) 
and  computed  wind  speeds  for  all  available  months.  The  random  error 
(approximated  by  root-mean-square-error,  RMSE)  and  bias  are  sufficiently 
small  that  they  are  not  expected  to  significantly  affect  the  statistical 
characteristics  of  the  hindcast  waves. 

8.  The  specific  comparisons  provide  descriptions  of  the  response 
of  the  wind  hindcast  procedure  (model)  to  the  supplied  input  (pressure 
gradients,  etc.)  for  specified  times  and  locations.  Since  the  model  can 
hindcast  winds  similar  to  measured  winds  in  the  specified  areas,  the 
model  is  expected  to  provide  reliable  results  in  areas  where  no  compara¬ 
tive  data  are  available. 


*  For  these  comparisons,  the  WIS  winds  were  transferred  to  the  10-m 
level  of  the  buoy  anemometers . 
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March  1973  at 
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DAYS  IN  APRIL  1973 

Figure  B14.  Comparison  of  NDBO  and  WIS  wind  data  for  April  1973  at  EB01 
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Figure  B15.  Comparison  of  NDBO  and  WIS  wind  data  for  August  1973  at 
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Figure  B16.  Comparison  of  NDBO  and  WIS  wind  data  for  September  1973  at  EB01 
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Figure  B19.  Comparison  of  NDBO  and'WIS  wind  data  for  December  1973  at  EB01 
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(  In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 

22  July  1977 ,  Subject:  Facoimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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